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CAN ONE MAKE A LASER OUT OF CARDBOARD? 


KRZYSZTOF BURDZY AND TVRTKO TADIC 


Abstract. We consider two dimensional and three dimensional semi-infinite tubes 
made of “Lambertian” material, so that the distribution of the direction of a reflected 
light ray has the density proportional to the cosine of the angle with the normal 
vector. If the light source is far away from the opening of the tube then the exiting 
rays are (approximately) collimated in two dimensions but are not collimated in three 
dimensions. An observer looking into the three dimensional tube will see “infinitely 
bright” spot at the center of vision. In other words, in three dimensions, the light 
brightness grows to infinity near the center as the light source moves away. 


1. Introduction 

We will examine the behavior of light rays in semi-inhnite tubes. The “cardboard” 
in the title of the paper refers to a material reflecting light according to the Lambertian 
distribution, to be described later in the introduction. The Lambertian distribution 
arises as the only physically possible reflection process in which reflected rays have 
random directions independent of the incidence angle (this follows from formula (2.3) 
in |ABS13] ). The “laser” effect refers to a possible collimation of light rays exiting the 
tube. We will show that if light rays are released far from the end of the tube and they 
reflect according to the Lambertian distribution then the exiting rays are collimated 
in two dimensions but are not collimated in three dimensions. So the answer to the 
question posed in the title is positive only in two dimensions. 

The three dimensional model does involve a singularity but of a milder type. We 
will show that an observer looking into the tube will see “inhnitely bright” spot at the 
center of vision. In other words, the light brightness grows to inhnity near the center 
as the light source moves away. 

The present project is a prelude to the study of Lambertian reflections in fractal 
domains. Some fractal domains have narrow channels and one would like to know 
how light travels within such channels. This article analyzes a toy model for the light 
behavior in a long thin channel. In future articles, we plan to extend this direction of 
research to light reflections in thorns with smooth boundaries and, ultimately, thorns 
with fractal boundaries. 

Our project is inspired by and related to a number of other projects. Lapidus and 
Niemeyer ( [LNlOl ILN13al ILN13b] ) considered billiards with the specular (classical) 
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reflection in fractal billiards. Comets et al. f |CPSV09l ICPSVlOal ICPSVlOb] ! studied 
random Lambertian reflections in smooth domains with irregular shapes. Angel et al. 
( |ABS13j ) showed that Lambertian reflectors could be approximated by deterministic 
reflectors. Evans f |Eva01j l studied a model of stochastic billiards were the reflection 
angle was uniform. 

We will describe the asymptotic behavior (angle and position) of the light ray when 
it reaches the end of the tube when the light source is far away. The motion of light 
rays along the tube is governed by a random walk. In order to hnd the exit position and 
angle of the light ray we need to hnd estimates for the distributions of undershoot and 
overshoot of a symmetric random walk. We will derive a number of explicit formulas 
using the Wiener-Hopf equation and various results from |Asm98[ IChowSG) IDon80( 
IEri70( IMik99[ |Rog71[ |Spi57| . See the book by Kyprianou |Kyp06| for an introduction 
to the topic. 

An intriguing and challenging aspect of the two dimensional model is that it leads 
to the “critical” case of the Central Limit Theorem. The model is associated with a 
random walk with steps that do not have a hnite variance but nevertheless the CLT 
holds (although we will not use this fact in our paper). 

The rest of the paper is organized as follows. We will present a more detailed overview 
of our main results in the next section. Section [3] contains a review of known results 
on random walks, Wiener-Hopf equation and related topics. We will derive there some 
new results needed later in the paper. Section 0] is devoted to the analysis of the 
two-dimensional model and hnally Section [5] presents results on the three dimensional 
tube. 


2. The model and main results 

We start with the description of Lambertian reflections of light. A physical surface 
is Lambertian if its apparent brightness does not depend on the angle at which the 
observer is looking at it. The Moon, in its full moon phase, is approximately Lambertian 
because it appears to be a globally flat surface to terrestrial observers despite being 
round. Lambertian reflections are also known as the Knudsen Law in the theory of 
gases. We will present the two-dimensional model in this section. See Section [5] for the 
three-dimensional case. 

Consider a set T) C with a smooth boundary. Suppose a light ray hits a point 
X G dD and reflects. The outgoing light ray travels at an angle 0 with the inward 
normal vector at x. The direction of the outgoing light ray is independent of the 
direction of the incoming light ray. The density of 0 is given by (see Figured]), 

^2 1 ) for 0 G (-7r/2,7r/2), 

|0, otherwise. 

The hrst part of the paper will be devoted to reflections in a semi-inhnite strip D = 
{{x,y) G : x < 0, 0 < i/ < 1}. We will assume that the light ray starts at (—s,0) 
for some s > 0 and travels in a direction which forms a random angle with the normal 
vector, with the density given by fl2.ip . The horizontal coordinate s of the starting 
point will play the role of the main parameter in our model. Whenever the light ray 
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Figure 1. Random reflection angle 0 and its density. 

hits the bonndary of Zi), it reflects according to the Lambertian scheme (see Fignre[2]). 
In particnlar, all reflection angles are jointly independent. At a certain time, the light 
ray will exit the strip throngh its opening dj-D := {{x,y) G : x = 0, 0 < y < 1}. 
Let ( 0 , 1 ( 5 ) be the exit point and let be the exit angle (see Fignre [2]). Onr main 
resnlt is concerned with the asymptotic behavior of the joint distribntion of (A^, Y() as 

s — )■ CXO. 



Figure 2. Starting point (—s,0), exit angle A^ and exit location ( 0 , 1 ( 5 ). 

The x-coordinates of the points of reflection constitnte a random walk. A step of 
this random walk is a symmetric random variable satisfying F{Xk > x) ~ cx~‘^ for 
X —)■ 00. An essential part of onr analysis is devoted to “nndershoot” and “overshoot”, 
defined informally as follows. The nndershoot Us is the horizontal distance from the 
last reflection point to drD. The overshoot Os is the difference between the size of the 
random walk step that goes beyond 0 and Us (rigorous definitions will be given below). 
One of our main results is the following simplified version of Theorem 14.101 

lim P ( \ < n I = tu^ 12, for t,u & [0,1]. 

s^oo yy logs Us + Os j 

Let 'll [a, h\ denote the uniform distribution on [a, h]. Our basic result on the limiting 
distribution for the exit angle A^ and exit location Ys, Theorem 14.131 says that, when 

s —)■ CXO, 

(a„K5)4(0,ii[0,i]). 

We use the results on overshoot and undershoot of the random walk to obtain more 
accurate information on the joint distribution of A^ and Ys in Theorem 14.141 For 
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t,ue [0,1], 

lim P 

S^OD 


'logcot lA^ 


log s 


< u, I As < 0,Ys < t 


u{l-{l-tf)/2, J = l, 
MtV2, I = -1. 


At this point we can answer the question posed in the title of the paper. We place 
the eye of the observer at approximately (0,?/) (see Figure [3]). The distribution of the 



Figure 3. Light rays arriving at the eye placed at (0,|/) = (0,2/3). 


light rays arriving at the eye is expressed in terms of A^ and given in Corollary 14.151 as 
follows. For u,y E [0,1] and e G (0, y) we have 


lim P 

S^OO 


'log cot I As 


logs 


< u, I As <0\Ys e {y 



u{l - y + e/2), 1 = 1, 

u{y-e/2), I =-1. 


The distribution is illustrated in Figure 01 Note the asymmetric singularity at 0. We 
continue the discussion of the two-dimensional results in Section 14.31 



Figure 4. Approximate distribution of As given Ys = y, with s = 1, 000 
and y = 2/3. 

We will discuss the three dimensional case in Section |5l The fundamental difference 
between two and three dimensional cases is that the asymptotic distribution of the 
direction of the light ray exiting the tube at a specihc point is degenerate in the two 
dimensional case and non-degenerate in the three dimensional case. We do not have 
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an explicit formula for the asymptotic exit distribution in the three dimensional case 
but we have some estimates. In three dimensions we have the following theorem of 
different nature. Let v* be the unit vector representing the direction of the light ray at 
the exit time assuming it leaves the tube at the center of the opening (see Section [5] for 
the rigorous dehnitions). Let B(r) = {[x,y,z) ■. + y'^ + = l,y'^ + < r'^,x > 0} 

denote a ball on the unit sphere. A somewhat informal statement of Theorem 15.101 is 


Theorem 2.1. For any 0 < ri < r 2 < 1, 


r’2 - ri 
27r‘^ 


Consider an observer at the center of the opening of the tube, looking towards the 
interior of the tube. The theorem says that small annuli at the center of the held of 
vision, with the area of magnitude 1/s^, receive about 1/s units of light. Hence, the 
apparent brightness is about s at the distance 1/s from the center, if the light source 
is s units away from the opening. This means that the surface of the tube does not 
appear to be Lambertian, i.e., the surface does not have uniform apparent brightness. 
This can be explained by the fact that not all parts of the surface of the tube receive 
the same amount of light. 


3. Review of stopped random walks 

In this section we establish notation that will be used throughout the paper, give some 
rigorous dehnitions, recall some known results and derive some theorems on general 
random walks, not necessarily those arising in the random rehection model. 

We will study a random walk {Sn,n > 0}, with S'o = 0 and S'„ = Sn-i + Xn for 
n > 1, where {Xn,n > 1} is an i.i.d. sequence. We will always assume that X„’s are 
continuous random variables. Some of the results stated in this paper might not be 
true for lattice variables. 

3.1. Renewal measures and ladder processes. The ascending ladder epochs {T^ : 
/c > 0} are dehned as 

(3.1) r+ = 0, 

r? = inf > Uli : Si > I , n > 1 . 

It is easy to see that {(T^ — T^j^,5'^+ — Srp+ ) : A; > 1} is an i.i.d. sequence. Let 
= Sj,+ — Srp+ ^ ioT k > 1. For n = 1,2,... we call R/" = Srp+ = 
ascending ladder heights. 

Similarly, we dehne the descending ladder epochs {TfT : A: > 0} by setting Tq” = 0 and 
T~ = inf{A; > T~_^ : Sk < Sj,- } for n > 1. The sequence {(Tr — Tr ,, Sj,- — Sj,- ) : 

A; > 1} is i.i.d. We let = Sj.- — Sj.- ^ for A; > 1 and call H~ = Sj,- = 
n > 1, the descending ladder heights. 

The following result can be found in |Don80] (see relations (4a) and (4b)). A more 
general sufficient and necessary condition for the hniteness of ladder step moments was 
given in |Chow86] . 
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Lemma 3.1. Suppose that E(Xi) = 0. 

(a) < oo then E[|Zj*"|] < oo. 

(b) E[X^] < oo if and only if E[Zf]E[—Zf] < oo. Moreover, 

E[X2] = 2E[Z+]E[-E-]. 


This immediately implies the following corollary. 

Corollary 3.2. If Xi is a symmetric random variable then E[Xf] < oo if and only if 
E[Z+] = E[-Zf] < oo. Moreover, 2E[Z+]‘^ = E[Xf]. 

We define renewal measures by 

OO 

u±(dx) = e dx). 

k=l 

One can show that for a measurable set A C M (see |Asm981 (2.4)]), 




This formula can be written in the following way. For a Borel set A C M, 

OO 

(3.2) U-(A) = P(Fo < 0,..., < 0, < 0, e A). 

n=0 


= E 




k=0 


Lemma 3.3. ITe have 


(3.3) 


P(Z+>«)=/ V(Xi> t-s)V-ids). 


Proof. Using the definition of Hf = Zf and fl3.2p . we obtain 
P(Z+ >t)= P(^1 >t)+ P(^1 < 0, ^2 > t) + P(^l < 0, ^2 < 0, ^3 > t) + . . . 

= [ P(Xi > t — s)x 

J —OO 

X (¥{So e ds) + P(Fo < 0, e ds) + P(S'o < 0, < 0, S 2 e ds) + .. 

= [ F{Xi>t-s)l]-{ds). 


□ 


The following result is the well known renewal theorem (see [DurlOl Sect. 3.4]; see 
|Eri70] for extensions). 

Theorem 3.4. For E[Z^] G (0, cx)] and all h > 0, 

ymUyiO.(])/« = ^ and Um 0+(|«, i + ft]) = 
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This implies that if G (0, cxo] then for all h > 0, 

(3.4) supI[J'''([f, f + h]) < cxo. 


t>o 


Definition 3.5. (a) For a function r : [0, oo) —)■ [0, oo), let 

CXD 

h sup{r(a;) : x e [{k — l)h, kh)}, 

k=l 

OO 

l’^{r) = h^inf{r(x) : x e [{k — l)h, kh)}. 


k=l 


We say that r{x) is directly Riemann integrable (d.R.i.) if lini/i^o= lim^^o 
and the limits are hnite. 

(b) Recall that the variation V{r) of r over [0, cxd) is dehned as 


sup E \r{xk) - r{xk 


-l)\y 


k=l 


where the supremum is taken over all sequences 0 < Xq < Xi < ■ ■ ■ < Xn- 
Remark 3.6. (i) It is elementary to check that 

I+{r)-I^{r) < V{r)h. 

ifT r{x)dx < oo, in the sense of the Lebesgue integral, then it is easy to see that 

poo poo 

I+{r) < / r{x)dx + — l'^{r)) < / r{x)dx + V{r)h, 

Jo Jo 

poo poo 

I-ir) > / r{x)dx — {I^{r) — I^{r)) > / r{x)dx— V{r)h. 

Jo Jo 

This implies that if V{r) < oo and J^r{x)dx < oo then r{x) is d.R.i. 

(ii) If r : [0, oo) —)■ [0, oo) is decreasing then it has a bounded variation. Hence, if r 
is decreasing and r(x)dx < oo then r(x) is d.R.i. 

(iii) Every d.R.i. function is necessarily bounded. Otherwise we would have /^(r) = 
oo for all h > 0. 

Lemma 3.7. Suppose that r : [0, oo) [0, oo) is d.R.i. and¥\Z^] G (0, oo]. Then 


(3.5) 


sup / r(s + a;)I[J'''((ia;) < oo, 

s>0 Jo 


(3.6) 


lim / r(s — x)I[J’''(dx) 

s^oo Jq 




r{s)ds. 


Proof. The claim fl3.6p can be found in [Purl01 (4.9)] or |Eri70l Thm. 3]). 

For fl3.5p we £x h G (0,1) and let Mi G (0, oo) be an upper bound for r (see Remark 
13.61 (hi)). By fl3.4l) . there exists M 2 > 0 such that I[J”'‘([t,f + h]) < M 2 for t > 0. Let 
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riQ = \_s/h\ + 1. Note that 0 < hqU — s < h. We have 


pnoh—s rinQi-Ki-ijn—s 

r{s + x)iJ~^{dx) = / r(s + x) I[J^((ia:) + ^ / r{s + x)V^{dx) 

Jo i._n J (nn+k)h—s 


/>{nQ+k+l)h—s 


k=0 '' [riQ+k)h-s 


< MiU''"(0, n^h — s) 

OO 

+ '^2 f sup{r(x) : X G [(no + k)h, {no + k + l)h)} 


k=0 

X U'''([(no + k)h — s, (no + k + l)h — s]) j 
<MiU+([0,l]) + (/^(r)/h)M 2 . 

The right hand side is finite and does not depend on s so fl3.5l) is true. 
For s > 0 we let 


□ 


(3.7) Ng = inf{n > 0 : S'„ > s}, 

Og = SN^-s, Ug = s-SN^_,. 

We call Og the overshoot and Ug the undershoot of the random walk Sn at s. We will 
also use the overshoot and undershoot of the ladder height process, defined by 

N+ = inf{n >0:H+ > s}, 

0+ = H+ - s, U+ = s- H+ ^. 

It is easy to see that 

(3.8) Og = 0+ and U+ < Ug. 

Lemma 3.8. //E[Zj^] = oo, then 0+ —)■ cxo and oo in probability as s ^ oo. 

Proof. We have 

OO 

P(f/^ <m) = ''^2 > s, s — m < < s) 

k=l 

= f ¥{Z+> s-y)V+{dy) <V+{[s-m,s]). 

J s—m 

The right hand side converges to 0 by Theorem 13.41 so Uf' —?■ oo in probability as 
s —)■ oo. 

A similar calculation yields 

OO 

(3.9) IP(0^ < m) = '^2 ^('S < ^k < s + m, < s) 

k=l 

P(s — y < Zf < s — y + m)l]^{dy). 

Note that 

P(s < Zf < s + m) < F{km < Zf < {k + l)m) + P((fc + l)m < Zf < {k + 2)m) 
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for km < s < {k + l)m. It follows that 



CX) « 

+ m)ds = / 

k=0 


p(k+l)m 
km 


P(s < < s + m)ds 


< ^ / \F{km < < {k + l)m) + P((/c + l)m < Z^ < {k + 2)m)jds 

k=0 


= m^P(/cm < Z^ < {k + l)m) + P((fc + l)m < Z^ < {k + 2)m)j < 2m < oo. 

k=0 


In other words, the function s —)■ P(s < Z^ < s + m) is integrable over [0, cxo). Since 

P(s < Z^ < s + m) = F{Z^ < s + m) — P(s > Z^), 

the function s i—)■ P(s < Z^ < s + m) is the difference of two monotone and bounded 
functions. It follows that this function has bounded variation. Since it is also integrable, 
it is d.R.i., by Remark [3]6] (i). Hence, by fl3.6p . 

lim [ P(s — y < Zf < s — y + m)iJ~^{dy) = 0. 

s^ooJq 

This and fl3.9p imply that 0+ —?■ oo in probability as s —)■ oo. □ 


Lemma 3.9. Make one of the following assumptions. 

(i) E[Xi] = 0, E[(Xf )2] < oo and E[(X+)2] = oo. 

(a) Xi is symmetric and E(X^) = 2E[(Xf )^] = oo. 

Then Of —)■ oo and Uf oo in probability as s ^ oo. 

Proof. Using Lemma [3.11 for case (i), or Corollary 13.21 for (ii) we obtain E[Z^] = oo. 
The claim follows from fl3.8p and Lemma 13.81 □ 

For functions /, g : (0, oo) —>■ M, we will write / g if \im^^^g{x)/f{x) = 1. 

Definition 3.10. For a function h : M’*' M’*' we say that it is regularly varying with 

exponent (index) a if 

(3.10) lim h{xt)/h{t) = a;“ 

t—^OO 

for a: > 0. A function h is called slowly varying if a = 0. 


Recall that h is a regularly varying function with index a if and only if it is of the 
form h{x) = x^L{x) where L is a slowly varying function. 

The following two results can be found in |Eri701 Thms. 6 and 7]. 


Theorem 3.11. Suppose that ¥{Hf > t) = t ^L{t), where L{t) is a slowly varying 
function. Then Of /s —)■ 0 and Uf /s —)■ 0 in probability as s ^ oo. 


Theorem 3.12. Suppose that Z has the distribution R[0,1]. If t ^ F{Hf > t) is 
regularly varying with index —1 and m{t) = j^F{Hf > x) dx then 


f m{Of) m{Uf) \ 
\ m(t) ’ m(t) J 


A{Z,Z), 


as t ^ oo. 
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The following theorem can be found in |Mik991 Thm. 1.2.4]. 


Theorem 3.13. Suppose that h : —)■ M"*" is regularly varying with index a < 0. 

Then for every a > 0, the limit in fld.lOh is uniform in x E [a, cxo). 

The following result, known as Potter’s Theorem, can be found in |BGT871 Thm. 1.5.6]. 


Theorem 3.14. Suppose that h : —)■ is regularly varying with index a. 

for any chosen <5 > 0 and e > 0 there exists to > 0 such that 


Ky) 

h{x) 


< 6 max 






for all t > to- 


Then 


Definition 3.15. A random variable Xi is called relatively stable if there exists a 
sequence of numbers a„ > 0, n = 0,1,..., such that Sn/an —?■ 1 in probability as 
n —)> oo. 


The following Relative Stability Theorem from |Rog71 Thm. 2] (see also Sect. 8.8 
in |BGT87] . especially Thm. 8.8.1) provides various characterizations of stable distri¬ 
butions. 


Theorem 3.16. If F{Xi > 0) = 1 then the following claims are eguivalent. 

(a) Ox/x —)■ 0 m probability as x ^ oo; 

(b) Ux/x ^ in probability as x ^ oo; 

(c) /;P(Xi > y)dy r'U L{x) where L{x) is a slowly varying function; 

(d) U’''([0, x)) = < x) ^ x/L{x) where L{x) is the same function as in (c); 

(e) Xi is relatively stable. 


The following result is taken from |Rog71 


Thm. 9]. 


Theorem 3.17. Suppose that Sn/an converge in distribution to a stable law with index 
a <2 for some seguence a„. If aF{Xi > 0) = 1 then Zf is relatively stable. 


Theorems 13.161 and 13.171 give the following result. 


Corollary 3.18. If Xi is symmetric and Sn/an converges to a normal distribution for 
some seguence an then 

(a) Zf and Zf are relatively stable; 

(b) U’''([0, x)) = U“((—x,0]) ~ x/L{x) where L{x) is a slowly varying function; 

(c) £F{Zf > y)dy = f^F(-Zf > y)dy r'U L{x) where L{x) is the same function as 
in (b). 

A sufficient condition for the convergence of Sn/an to a normal distribution is con¬ 
tained in the following very general theorem (see |Mik991 Gor. 1.4.8]). 

Theorem 3.19. Let P(|l^i| > t) = t~‘^L{t) where L is a slowly varying function. Then 
Xi is in the domain of the attraction of the normal distribution. 


Lemma 3.20. Suppose that Xi is a continuous symmetric random variable and t i-A 
P(Xi > t) is regularly varying with index —2. The following claims hold: 























CAN ONE MAKE A LASER OUT OF CARDBOARD? 


11 


(a) IE[Xf] = oo and = oo; 

(b) U+([0,t)) ~ t/L{t) where L{t) = f*F(Z^ > x) dx is slowly varying. 

(c) l[J’''(t(is)/I[J+([0,t)) converges weakly to the Lebesgue measure on [0,^4] for any A > 

0 . 

(d) The following limit is uniform in s > 0, 

F{X,>t{l + s)) _ 1 

P(Xi>t) (1 + S)2- 

Moreover, for every £ > 0 there exists to > 0 such that 


(3.11) 


P(Xi>t(l + g)) ^ 1 

F{Xi>t) -(l + s)2-^’ 


for all t > to- 


Proof, (a) It is easy to check that if t i-A P(Xi > t) is regularly varying with index —2 
then E[Xf] = oo. Part (a) then follows from Corollary 13.21 

(b) By Theorem 13.191 the assumptions of Corollary 13.181 are satished. We can take 
the same slowly varying function function L{t) = J^F(Zf > x) dx in Theorem 13.161 fcl 
and Corollary 13.181 fexcept in this case the hrst positive step is denoted Zf). Part (b) 
of the lemma now follows from Corollary 13.181 (b). 

(c) It follows from Corollary 13.181 (b) that 

U’''(t[0,a]) at/L{at) 

U+([0,t)) ^ t/L{t) ^ 

for all a G [0, A]. This implies part (c) of the lemma. 

(d) The claims follow from Theorems 13.131 and 13.141 □ 


Lemma 3.21. If Xi is a continuous symmetric random variable and t i-A P(Xi > t) is 
regularly varying with index —2 then 

(3.12) - 


> t)u+([0,f)) 


= 1 . 


Proof. From formula fl3.3p we have 

poo 

P(^ 1 + >t)= P(Xi > t + s)U+(ds) 


= P(Xi >f)U+([0,f)) 

The change of variable s = tu gives us 

Hzt >t) 


r°°P(Xi>f + s) U+{ds) 

0 F{X,>t) U+([0,f))' 

P(Xi > f(l + «)) lJ+{tdu) 


P(Xi>f)U+([0,f)) Jo F{X,>t) U+([0,f))- 

It will suffice to show that 


lim 

t—>00 


P(Xi > t(l +«)) \J+{tdu) 


F{X,>t) U+([0,f)) Jo (1 + s)' 


ds = 1. 
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It is not hard to see, using Lemma 13.201 (c) and (d) that 

P(Xi > t{l + m)) U+(t du) 


(3.13) 


lim 

t—^OO 


P(Xi>t) u+([0,t)) 




ds. 


According to (I3.1ip . we can pick to > 0 such that P(Xi > t(l + s))/P(Xi > t) < 
H{s) := (1 + for t > to- For ^ > to, using the integration by parts formula, 


P(Xi > t(l + m)) m{tdu) 
P(Xi > t) U+([0,t)) 


< 


H{u 


i]~^{tdu) 

fFoTO) 


(3.14) 


= H(u) 


D+(|0, 111)) 


U+(|0,«)) 


2=A 


f°° w(|o,fa)) 

L D+([0,i)) 


H'{u) du. 


Recall the slowly varying function function L(t) = JqP(Xi > x)dx and use Lemma 
13.201 (b) to find to such that 


‘ <uyio,())< 


2L(t) ^ ^ ^ L{ty 

for t > to. From the fact that L(t) is increasing for t > A we obtain 

, < 4., 

D+([0,4)) “ L{ts) - 

Since H{s) ~ and H'{s) ~ (I3.14p is bounded by 

poo 

AAH{A)- / AsH'{s)du. 

J A 

Choose an arbitrarily small £ > 0 and pick A large enough so that 


P(Xi > t(l +m)) \]+{tdu) 
P(Xi > t) U+([0,t)) 


< £ 


and 


< A 


(1 + s)' 


ds < e. 


It follows from this and fl3.13p that for large t, 

>t{l + u)) V+{tdu) 


- 1 


< 


P(Xi>t) u+([0,t)) 
p(Xi > t(i + m)) m{tdu) 

) P(Xi > t) U+([0,t)) 


(1 + .)^ 


■ ds 


Ja (1 + 

The claim now follows. 


ds T 


P(Xi > t(l +m)) l]+{tdu) 
P(Xi > t) U+([0,t)) 


E -\~ S S. 


□ 


Corollary 3.22. Suppose that Xi is a continuous symmetric random variable and 
t I—)■ P(Xi > t) is regularly varying with index —2. Then t e-?> P(Z^ > t) is regularly 
varying with index —1. 


Proof. The corollary follows from Lemmas 13.201 (b) and 13.211 


□ 
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Lemma 3.23. Assume that Xi is a continuous symmetric random variable and t i—)■ 
P(Xi > t) is regularly varying with index —2. Then 


(3.15) 


m 


(t) := / > x) dx 



2xP{Xi > x) dx. 


Proof. Using Lemmas 13.201 fb) and 13.211 we get 


(3.16) 


lim 


P(Z+ > t) jl P(Z+ > x) dx 


= 1 . 


tP{Xi > t) 

Since = oo, Corollary 13.21 yields E[Zf] = oo. Therefore, m{t) f oo and 

n{t) := fgxP(Xi > x) dx f oo as t —>■ oo. Hence, we can use rHopital’s rule to 

m'^(t)/n{t) , and we get 

> x) 2m’{t)m{t) 


calculate the limit lim 


lim 


= lim 


/JxP(Xi > x) dx n'{t) 


= 2 . 


The last equality follows from (I3.16p . This easily implies the lemma. 


□ 


The lemma easily implies the following corollary. 


Corollary 3.24. Suppose that Xi is a continuous symmetric random variable such that 
limi^oo > t) = c E (0,cx)). Then 

F{Zf > f) ~ f > x) dx ~ V^-^/logf. 

t^J\ogt Jo 

Lemma 3.25. Suppose that Xi is a symmetric random variable such that t i-A P{Xi > 
t) is regularly varying with index —2. 

(a) Os/s —0 in probability when s ^ oo. 

(b) If m{t) = \JJq 2xP(Xi > x) dx then m{Ot)/m{t) ^[0,1] in distribution when 

t oo. 


Proof, (a) Note that Hf = Zf. Part (a) follows from Corollarv l3.221 fl3.8p and Theorem 

EZH 

(b) Once again, we will use the fact that Hf = Zf. Part (b) follows from fl3.8l) . 
fl3.15l) , Corollary 13.221 and Theorem 13.121 □ 

Lemma 3.26. If Xi is a symmetric random variable such that linit^oo > t) = 

c e (0, oo) then y^log Ot/ logf — )■ lf[0, 1] in distribution as t ^ oo. 

Proof. The lemma follows from Corollary 13.241 and Lemma [3.251 (b). □ 


3.2. Wiener-Hopf equation. The Wiener-Hopf integral equation is 
(3.17) PU(s) = ( 7 (s)+ r W{s-y)F{dy), 

J —OO 

where W : [0, oo) —)■ M is an unknown function. The function g : [0, oo) —?■ M and the 
probability distribution P on M are given. 

We will make the following assumptions, common in this context. 
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• g{x) > 0 for all X > 0 and sup 2 .g[Q^„] g{x) < oo for all a > 0 . 

• F is a probability measure with a well defined mean. 

• We will consider only positive solutions to fl3.17p . i.e., W{x) >0 for all x > 0. 

If = 0, then we call the equation homogeneous. Spitzer has shown in |Spi57| that, 
in general, there is no uniqueness for solutions to the homogeneous equation. However, 
uniqueness holds if F is concentrated on [0, oo); see |Durl0] . 

In this paper, F in fl3.17p will be the distribution of Xi. For s > 0 we dehne 

Mk = maxIFj : 0 < j < n}, 

OO 

(3.18) fFmin(s) = ^E(5f(s - S'fc)l(M;,<s))- 

k=0 

We will need the following result from |Asm981 Cor. 3.1]. 

Theorem 3.27. Any solution W of the equation fl3.17p is of the form W = Wmin + hFo; 
where Wq is a solution to the homogeneous equation. The function IFmin(s) defined in 
fl3.18l) is the minimal solution. 

Once again we quote a result from |Asm98l Prop. 3.3]. 

Theorem 3.28. For IFmin we have 

fPmin(s)= [ G{s-x)U^{dx), 

Jo 

where G{s) = J°^g{s - y)V-{dy). 


Lemma 3.29. Let F be the probability distribution function of a symmetric random 
variable such that 


\x\F{dx) < oo. 


\x\^F{dx) = oo. 


and assume that for all s > 0 , 

g{s) < d{s)r{s), 

where d and r are directly Riemann integrable and d is non-increasing. Then, Wmin the 
minimal solution to the equation (I3.17p has the property that 

lim HAiin(s) = 0 . 

S—^OO 

Proof. Let {Xk} be a sequence of i.i.d. random variables with distribution F. Since Xi 
is a symmetric random variable, I[J’'‘((ix) = U“(—dx) and we set U = U’'". Using the 
notation of Theorem 13.281 

/ O poo 

g{s - y)\J-{dy) < d{s) / r{s + y)V{dy) < Gd{s), 

■OO J 0 

where G := supg>o/o°° + y)^{dy). It follows from our assumptions and Corollary 

13.21 that ¥,Zf = oo. Hence, we can apply fl3.5p to see that the constant G is finite. By 
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fl3.6p and Theorem I3.28[ 

0<hh^in(s)= rG{s 
Jo 


x)V{dx)<C [ d{s — x)iJ{dx) ^ 0, 
Jo 


when s —?■ oo. 


□ 


Corollary 3.30. Suppose that F is the probability distribution function of a symmetric 
random variable and 1 — F{x) is regularly varying with index —2. Assume that there 
exist C > 0 and a > 0 such that 

( 319 ) 9 ( 0 ) ^ 

Then lirng^oo hhmin(5) = 0, where hhmm is the minimal solution to the eguation fl3.17p . 


Proof. Let d{s) 


= 1+^/2 • Then for all s > 0, 

^ 1 + 52+a ^ (^l + sl+«/2 j 


d{s)r{s). 


Since s i-A is a decreasing and Lebesgne integrable fnnction on [0,oo), it is 

directly Riemann integrable, by Remark 13.61 (ii). The claim now follows from Lemma 


[3291 


□ 


Remark 3.31. Corollary I3.30l mav not hold for a = 0, as the following example shows. 
Let F be the cnmnlative distribntion fnnction of a symmetric random variable with 
1 — F{x) = for X > 0. Let Ng denote the stopping time dehned in (13.7p for the 
random walk with the step distribntion F. In this case we have Sn/n —?• 0 a.s., and 
by the Chnng-Fnchs Theorem the random walk is recnrrent, hence Ng < 00 , a.s. By 
Theorem 13.271 the eqnation fl3.17p with g{s) = 1 — F{s) for s > 0 has the minimal 
solntion 

OO OO 

Wmin{s) = 1 = P(iV^ < 00 ) = ^ F{Ng = n) = ^ F{Xk > s- Sk-i, Mk-i < s) 

n=l k=l 

00 

k=l 


4. Two-dimensional model 

Recall the two dimensional model from Section |2l First we will review some proper¬ 
ties of the random angle 0 with the density fnnction given by fl2.ip . The cnmnlative 
distribntion fnnction F© is eqnal to 

fo, t<vr/2, 

Fe{t) = P(0 < t) = < J{smt + 1), t e {-7rf2,7rf2), 

[l, t>TT/2. 

Note that F©^ : (0,1) ^ (—7r/2,7r/2) is given by 

Fe\y) = sin-^(2|/- 1). 
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If V has the distribution 'U(—1,1) then it is easy to check that the following equalities 
hold in the sense of distribution, 


(4.1) 


sin Q = V, cos 0 = Vl — V^, tan 0 = 

VI -142 


tan© b+tan©,l) 



Figure 5. Step of the random reflection. 


In the random reflection model described in Section |21 if the ray is reflected at the 
point (x, u) then its next reflection point will be at {x + tan 0,1 — u), where u G {0,1} 
and 0 has the density given by fl2.ip (see Figure |5]). 

Let {0n} be a sequence of i.i.d. random variables with density given by (12.1 p and 
set Xn = tan0„. We dehne a random walk by setting Fq = 0 and Sn = Sn-i + Xn for 
n > 1. Recall that Ng = inf{n > 0 : > s}. Then the trajectory of the light ray 

described in Section [2] consists of 

(i) line segments [(-s + Sn-i, (1 - (-l)”"V/2), {-s + Sn, (1 - (-l)”)/2)] for n < W; 

(ii) line segment between {—s+Sns-i, (1 —(—1)^“~^)/2) and {—s+Sm^, (1 —(—1)^®)/2)- 
In view of dn]), the representation (i)-(ii) given above can start alternatively with a 
sequence {K,} of i.i.d. 1X(—1,1) random variables and Xn = 14/■\/l — V^. 

Definition 4.1. We dehne to be the angle between the exiting ray given in (ii) and 
the inward normal to the right edge [(0,0), (0,1)]. We let Yg denote the ^/-coordinate 
of the point where the ray exits the tube through the right edge (see Figure |2]). 

Lemma 4.2. (a) The cumulative distribution function of Xi is Fx,{x) = I + 
and its density is fxi{x) = (1 -|- x'^)~^Ff2 for x G M. 

(b) E[A:„] = 0 and E[Xl] = oo. 

(c) The random walk {Rn} is neighborhood recurrent. 


Proof. Part (a) follows from an elementary calculation, 
(b) Since 


E|W„| = E 


|W 


X 


7-1 yr 


X 


= - dx = 
22 


X 


vT 


: dx = 1, 


x^ 


we must have EX„ = 0 by symmetry. The second moment is inhnite because 

2x2 


EK] = 


L^(i+x^rn 


dx = oo. 
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(c) Since EXi = 0, the strong law of large nnmbers shows that Sn/n —)■ 0, a.s. This 
also holds in probability so the Chnng-Fnchs Theorem for random walks implies that 
{Sn} is a neighborhood recurrent random walk (see [DnrlOl Thms. 4.2.1 and 4.2.7]). □ 

It follows from Lemma 1121(c) that the light ray will hit the line {x = 0}, a.s. In other 
words, with probability 1, the light ray will exit the tube (strip) through the right edge. 


Lemma 4.3. For s > 0 and ti > t 2 > 0 we have 


am l.mhALLi^4. 

V \ „/+0 +2 S^oo P(Xi > 5 /^ 2 ) 


^2 


P(Xi > s/h 

Moreover, for all t G (0,1] and all s > 0 we have 
(4.3) 0 < P(Xi > s/t) - f2p(Xi > s) < 


(1 + s2)2- 


Proof. Both claims in (I4.2jl follow easily from the following identity, 


P(Xi > s/F) 
P(Xi > s/t2) 


1 

s 

*1 

2 

2\/tf+s2 

^ytl+s^{^//q+s'2+s) 

1 

s 

^2 

2 

2\/*2+®^ 


''1 

\/tl + s'^ 

f^/tl + s^ + s) 

to to 

+ s2 

'{y/ff + s2 + s) 


We have 


P(Xi > s/t) - t¥(Xi > s) 
1 


= t 


(4.4) 


+ s) \/l + s2(\/l + s2 + s) 

1 + ^ 

_ ^2\ _ \/G+.s^+\/l+s^ _ 

(f2 + §2 s^/F + s2)(l + s2 + S\/l + s2) 


Clearly, the expression in (14.4p is non-negative for t G (0,1] and s > 0, f^(l — f^) < 1 
and since s > 0 we have 1 -|- < 2. For s > 1, 


1 


< 


(f2 4 - §2 _|_ Sy/W+~^){1 + §2 + s\/l -|- s2) (1 + s2)2’ 

SO by (Ol) we get P(Xi > s/t) - t^¥{Xi > s) < For s < 1 we have 

4 


P(Xi > s/t) - FF{Xi > s)< P(Xi > s/t) < 1 < 


(1 + s2)2' 


In either case, fl4.3p holds. 


□ 


Since E[Wf] = 00 , we conclude from (13.8p and Lemma [3.91 that, in probability, when 
s ^ 00 , 

(4.5) Os = Sn^ — s —)■ 00 and Us = s — —)■ 00 . 
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Let 

u(s, t)=¥ ^ > s - 5jv,-i) 

CX> 

= j2nxN. > {s-SM.-i)/t,N, = k) 

k=l 

oo 

(4.6) = > (s - Sk-i)/t,Mk-i < s). 

k=l 

Note that 

(4.7) u(s, 1) = 1, 

because u(s, 1) is the probability that the random walk {Sk} will take a value greater 
than s and by Lemma [4.21 (c), this probability is 1. 

Lemma 4.4. u(s, t) is the minimal solution to the Wiener-Hopf equation 
W{s,t) = P{X^> s/t)+ [ W{s-y,t)¥xAdy)dy. 

Proof. Fix t > 0 and let g{s) = P{Xi > s/t). Formula fl4.6p and the Markov property 
imply that u{s,t) = ~ Sk)lM^,<s)- By Theorem 13.271 the function u(s,t) is 

the minimal solution to the Wiener-Hopf equation. □ 


Lemma 4.5. The function 

u(s, t) := u(s, t) — f^u(s, 1) = u(s, t) — t^ 
is a solution to the following equation 


hF(s, t) = P(Xi > s/t) - t^P{Xr >s)+ [ W(s - y, t)fx{y) dy. 

J —OO 


Moreover, this is the minimal solution to this equation, that is, for every (positive) 
solution W of this equation we have u(s,t) < W{s,f). 

Proof. We have 

u(s, f) — t^u{s, 1) 


oo OO 

= ^ P(Xfc > (s - Sk-i)lt, Mk-i < s) - ^ t^^{Xk > s- Sk-i, Mk-i < s) 

k=l k=l 

oo 

= ^ P(Xfc > (s - Sk-i)lt, Mk-i < s) - t^P(Wfc > s - Sk-i, Mk-i < s) 

k=l 

= P(Xi > s/t, Mk -1 < s) — t^¥{Xi > s, Mk-i < s) 


+ ^ F{Xk > (s - Sk-i)/t, Mk-i <s)- f^P(Xfc > s - Sk-i, Mk-i < s). 

k=2 
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Setting g{s) := P(Xi > s/t) — t‘^F{Xi > s), we obtain, 

(4.8) P(Xfc > (s - Sk-i)/t, Mk-i <s)- i^P(Xfc > s - Sk-u Mk-i < s) 

= [ P(Xfc > (s - y)/t) - t‘^¥{Xk > s- y)¥{Sk-i G dy, Mk -2 < s) 

) 

g{s - i/)P(S'fe_i e dy, Mk -2 < s) 


Hence, u(s,t) = YlT=i^i9i^ ~ Sk-i)l{Mk-i<s)), and from (13.181) and Theorem 13.271 we 
know that this is the minimal solution to the Wiener-Hopf equation in the statement 
of the lemma. □ 


Lemma 4.6. For a fixed t E (0,1] we have 

lim u(s, t) = lim u(s, t) — f^u(s, 1) = 0. 

s—)-oo s—)-oo 


Proof. By (14.3p . P{Xi > s/t)—FP{Xi > s) < 4(l + s^)“^. Since 2(l + s^)“^ ds < cxo 
and s I—)■ 2(1 + s^)“^ is decreasing. Remark iLOl (ii) shows that this function is directly 
Riemann integrable. This implies that 4(l + s^)“^ is a product of two decreasing directly 
Riemann integrable functions. The lemma now follows from Lemmas 13.291 and 14.51 □ 


Lemma 4.7. (a) We have lims^oo u(s, t) = f^ for all t G [0,1]. 

(b) Suppose that Z has the distribution 1X(0,1). The following limits hold in distri¬ 
bution, as s ^ oo, 


(4.9) 

(4.10) 


f/.s 


0, + Us 




o. ^ 1 - ^fz' 


(c) The following limits hold in probability, 


lim 

S^OO 


log log O, 
lim -^- 

s^oo log S 

\og{Us + Os) - logO^ 
log s 


0 , 

0 . 


Proof, (a) Recall from fl4.7p that u(s, 1) = 1 and apply Lemma [4.61 

(b) Since is the cumulative distribution function of \pZ and, by dehnition, u(s, t) = 
(^ o^+u — ^)’ follows. The formula in fl4.10p follows easily from fl4.9p . 

(c) Take the logarithm of the left hand side in fl4.9p (resp. (I4.10p ) and divide by 

logs. The logarithm of the right hand side of each (14. 9 p and (I4.10p . divided by logs, 
converges to 0 in distribution. □ 
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Corollary 4.8. Suppose that¥{Xi > t) = ^ distribution U{0, 1). 

The following limits hold in distribution, as s ^ oo, 


(4.11) 


(4.12) 


I log Us \og Os 


log s ’ y log s 


l\og{Us + Os) \og Os 


log s 


log s 






Proof. The corollary follows from Lemma [3.261 and Corollary 14.71 (c) 


□ 


4.1. Asymptotic independence of exit characteristics. From the intuitive point 
of view, one would expect that when s, the distance from the light source to the right 
edge of the strip, is large then the following random variables would be approximately 
independent: the size of the undershoot, the ratio of the undershoot and overshoot, and 
the last side (upper of lower) visited by the light ray before the exit from the strip. We 
will prove that this is actually true. The idea of the proof is natural but its rigorous 
presentation requires extensive formulas. 


Lemma 4.9. For t G [0,1] and j = 0, 1, 


lim P 

S^OQ 


' log Us 
log s 


< t, 12 n{Ns) = j] = -t- 


Proof. We set Sn = and = Ylk= 2 ^k- We dehne Ns = infjn : Sn > s}, 

A' = inf{n : S!^ > s} and Us = s — U'^ = s — The dehnitions of Os and 

Og are analogous. We have 

'log Us 


p 


= p 


logs 
logU 


< t,Ns e 2N 


P 


logs 
log A 


log s 


< Ns G 2N, min{0„ O^, Us, U'}/2 > |Ai|^ 
< t\ Ns G 2N, min{0„ O^, U'g}/2 < |Xi 


On the event minjO^, O'g, Us, C'}/2 > |Ai| we have U^ = Us — Xi and = A' + 1. 
Hence we have 
'log A, 


P 


< P 


< t,Ns G 2N 


< P 


log s 

/ log-(!7: + X 

V logs 

p(min{o..o:,t;.,t;;}/2< |X,|) 
‘“sC <^2^iog(i + Xi/t;:) 


< e, a: G 2N - 1, min{0„ O^, A„ A:}/2 > | Ai 


log s log s 

+ P(min{0„0;,A„A;}/2< |Ai|) 


, a; G 2N - 1, min{0„ O', A„ U'g}/2 > | Ai|^ 
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< P 


\ogU'^, , log(3/2) 


<V + 


log s log s 

+ P(min{0„0:,t/„t/;}/2< |Xi|) 


, iV; e 2N - 1, min{0„ O^, 17:}/2 > \X^\^ 


< P 




log s log s 

+ P(min{0„0;,t/„t/;}/2< |Xi|). 

Since (S'„) and {S'^) have the same distribution, 


P 


1^ 

log s 


< t^iV' e 2N - 1 ) = P 


log^. 

log s 


< t^,Ns e 2N - 1 


Hence, subtracting from the previous inequality we get 


P 


log 17, 


log s 
< p I < 


< e 2N 


P 


log 17, 
log s 


< 7f, e 2N - 1 


2 . ^OgUi 


log s 


<^2^1og(3/2) 


log s 


, iV' e 2N - + P (min{0„ 0'„ f/„ U'}/2 < |Xi|) 


- ^ O' ^ + ^TogfO ^ ^ < 1^11) • 

It follows from Lemma 13.261 and (14.111) that min{0,. O', 17,, f/'} —)■ cxo in probability as 
s ^ cxo. The hrst term on the right hand side of the last formula goes to 0 as s —)■ cxo 
in view of fld.lip . Thus, 


lim sup P 

S—)-oo 


log 17, 
log s 


< t^,Ns e 2M 


P 


log 17, 


We can show in a similar manner that 


lim sup P 

S^OO 


log 17, 
log s 


< t^,Ns e 2M - 1 


log s 
log 17, 


< W e 2M- 1 = 0. 


P 


The claim now follows from the fact that 


lim P 

S^OO 


log 17, 
log s 


< iV, e 2N - 1 + P 


log s 

log Us 
log s 


< iV, e 2N = 0. 


< t^iV, e 2N =t 


□ 


The following is one of our main results. 
Theorem 4.10. For t, u G [0, 1], j = 0,1, 

17, 


lim P 

S^OO 


'log 17s 
logs 


<t, 


Us + 0, 


< V, l2n{N,) = j ] = -tv^ 


Proof. First note that 
P 


logf/s^^ Us ^ ^ ^ 

— 77— r~PV — 12n(-1Vs) — 0 


log s 


Us+Os 
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= ^ { Us <s\ 


2 Us 


= P s - Sns-1 < S , 


Us+Os 

£2 3- Sms-1 


< V, l2n{Ns) — 0 


X 


Ns 


< V, 12n(Xs) — 0 


( s - < /, ^ <V, Ns = 2k 


k=l 

oo 

J2p 

k=l 

PS 


X 


Ns 


3 — S2k-l 


P 


V 

3 — U 

V 


^ 3 — 3* < S2k-l < S, S2k-2 < S, . . . , S'! < S 


< Xi J P {S 2 k-i £ du, M 2 k-i < 3 ) 


k=l 


We have 


^:=P 


'^ogUs Us n n 

-j - < t, < V, 12n(^J = 0 

log s Us + Os 


— T^P 


log 3 


< t, 12h(-^s) — 0 


P j S U ^ ^ ^ ^2 


< Xi - U^P (s - M < Xi 


IP (*S' 2 fe-i G du, M 2 k-i < 3 ) 


k=l 


Lemma [4.31 implies that > 0. It follows from (14.Sh that 

OO 

^P(S'fc_i e du,Mk-i < 3 ) 

k=l 

= u(s,t), 

where u is dehned in Proposition 14.51 By Theorem 14.61 lim^^oo u(s, v) = 0. The 
theorem now follows from Lemma 14.91 □ 


41. < 


P 


3 — u 


<Xi 


— v 


*P(S-M < Xi) 


We record a few variants and corollaries of the last theorem. They follow easily from 
Lemma 14.71 (c) and Theorem 14.101 


Corollary 4.11. For f, n G [0, 1], j = 0,1, 


(4.13) 

(4.14) 


lim P 

S^OO 


lim ] 

S^OO 


'log Os 


<t, 


Us 


log S Us + Os 


< V, 12n(^.) = j] = 


'\og{Us + Os 


log 3 


< t 


Us 


' Us + Os 


< V, 12n{Ns) = i] = 


Us , , \ 1 2 

.- I tt—TT - I 2 N at, =j]=-v, 

s^oo \Us + Os / 2 


lim P 
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(4.15) lim P (iV, e 2N) = lim P (iV, e 2M - 1) 


1 

2 


4.2. Exit angle and position. We introduced the exit angle and position Yg in 
Definition 14.11 Now we will describe their joint distribution. Recall that (s + : 0 < 

n < Ng — 1) are x-coordinates of the reflection points of the light ray inside the tube 
before the exit time. 

Lemma 4.12. For s > 0 we have 



(4.16) 


Proof. If Ng is even then the last reflection happened on the upper boundary of the 
tube and the angle is negative. Hence, = — cot“^(Os + Ug). One can use similar 
triangles (see Figure |6]) to show that Yg = 



Figure 6. The case when Ng is even and A^ is negative. 


The case when Ng G 2N — 1 can be dealt with in a similar way. 


□ 


Theorem 4.13. (A^, W) (0,1I[0,1]) in distribution as s ^ oo. 

Proof. Recall from (14.51) that —)■ oo in probability. Therefore A^ —)■ 0 in 

probability as s —)■ oo. 

It remains to show that W —)■ R[0,1]. We use fl4.14p . fl4.15p and fl4.16p to see that 


P(5A <t)= P(W < t, A^. e 2N - 1) + P(w <t,Nge 2N) 



□ 
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Theorem 4.14. Fort,v G [0,1], 

< tjAs < o,y; < T I = 


lim P 

S—^OO 


'log cot I As 


log s 


tv^/2, if j = -1. 


Proof. It follows from Lemma [4.121 that log cot |As| = log(Os + Ug) and 

{Ng e2N-l,-^^ <v}, ifj = -l, 


{jAs < 0, Ts < w} = 


{Ng G 2N, 1 - 


Os-\-Us 

Us 

Os-\-Us 


<v}, ifj = l. 


The theorem now follows from fl4.13l) . 

Corollary 4.15. For t,y E [0,1] and £ G (0, y) we have 


□ 


lim P 

S—^OO 


'log cot I As 


log s 


< t,j Ag < 0 I Ts e (2/ - i/j = 


t{l-p + £/2), ifj = l, 
t{p-£f2), ifj = -l. 


Proof. Theorem 14.131 shows that lims^ooIP(^s £ {v ~ ^^v]) exists and is non-zero. This 
and Theorems 14.131 and 14.141 can be used to derive the asymptotic formula for the 
conditional probability. □ 


4.3. Discussion of the results. Theorem 14.131 savs that in the limit (i.e., when the 
light source is inhnitely far away), the light rays exit the two-dimensional tube hori¬ 
zontally, and they are equally likely to exit at any point of the right edge. 

Next we discuss the direction from which light rays arrive at an eye located at a point 
( 0 , 1 /) (see Figure 13]) • Corollary 14.151 savs that for large s, 

where V has the uniform distribution ^^[0,1] and R is an independent random variable 
with P(i? = 1) = y and P(i? = —1) = I — y. We can “solve for A*” to derive the 
following purely heuristic formula, 

A^ ^ Rcot~^ ■ 

Approximately y proportion of light arrives from the lower side (yellow rays in Figure 
13]), while the remaining rays arrive from the upper side of the tube (orange rays in 

Figure 13]) • The histogram in Figured represents a simulation of i?cot“^ 


5. Three-dimensional model 

This part of the paper will be devoted to light reflections within a three-dimensional 
semi-inhnite cylinder C = {{x,y,z) G + y'^ = l,x < Q} (see Figured). In this 

case, the exiting light rays are not asymptotically parallel when the light source moves 
to inhnity. So the three-dimensional model is less degenerate than the two-dimensional 
model. In this case, our results are less complete than those in the two-dimensional 
case. The reason is that deriving explicit formulas for this model is hard—this is a well 
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known difficulty with models related to the Wiener-Hopf equation (see Section 6.5 of 
|Kyp06| , and especially subsection 6.5.4). 



Figure 7. Cylinder 

We will assume that the light ray starts at s := (—s, 0, —1). At the initial time and 
whenever the light ray hits the boundary of C, it reflects according to the Lambertian 
scheme, i.e., 

(i) the outgoing light ray forms a random angle 0 with the normal to the tangent 
plane, 

(ii) 0 is a random variable with density fl2.ip . 

(iii) the projection of the outgoing ray onto the tangent plane forms a random angle 
$ with the line parallel to the x-axis (see Figure [8]), 

(iv) $ has the distribution If)—7r/2,7r/2] and is independent of 0. 

The consecutive reflection directions are jointly independent. 



Figure 8. Reflection with respect to the tangent plane 


Consider the light ray leaving the starting point {x, 0, —1). The tangent plane to the 
cylinder C at that point is {z = —1} and the ray starts moving along the line parallel 
to the vector 

(5.1) (i?sin0cos<h,i?sin0sin<h,i?cos0), R>0. 


Lemma 5.1. Given 0 and <h, the distance to the next reflection point is 


(5.2) 


R 


2cos0 

sin^ $ sin^ 0 + cos^ 0 


Proof. We need to hnd a point (i?sin0cos<h + x, Rsin 0 sin $, Rcos 0 
cylinder {ifl + = 1}. A straightforward calculation yields the formula. 


1) on the 
□ 
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Lemma 5.2. If the light ray reflects at the point (xo,yo,Zo) and 0 and $ are given 
then the next reflection point will occur at 

(/?sin 0 cos $ + xq , R{—zo sin 0 sin $ — yo cos 0), R{yosm 0 sin $ — zq cos 0)), 

where R is given by fl5.2p . 


Proof. Note that (xo,yo,Zo) = (xqjO, —1)^4 where A is the matrix representing the 
rotation operator about the x-axis and given by 


A = 


1 0 0 
0 -zo -yo 
0 yo -Zo 


In view of fl5.ip - fl5.2p . if the light ray starts at (xq, |/o, ^o); then the next reflection point 
will be at the point (i? sin 0 cos <h + xq, -R sin 0 sin <h, i? cos 0)^4. □ 


Next we establish notation for the process of reflection points inside the cylinder 
C. Recall that the light ray starts at s = (—s,0, —1). The reflection points will be 
{Sk + s}fc>o where {Sk}k>o is a random walk defined as follows. Let (0^, be an 

i.i.d. sequence such that 

• 0fc has the density |cos0 on [—7r/2,7r/2] for all k] 

• is distributed as 'U[—7r/2,7r/2] for all fc; 

• all random variables in the union of the families {Qk)k>o and (4)fc)fc>o are inde¬ 
pendent. 

Set So = (0, 0, —1), 

^ 2 cos 0fc 2 cos 0fc 

^ sin^ sin^ 0^ -|- cos^ 0^ 1 — cos^ sin^ 0^ ’ 

and define Sk = {Sf, Si, S^) for /c > 1 by 

(5.3) Sf = Sl_i + Rk cos ^k An Qk, 

(5.4) SI = Sl_, + Rk(-Sl_, sin sin 0^ - Sl_, cos 0^), 

(5.5) SI = Sl_^ + Rk(Sl_^ sin 4)^ sin Ok - Sl_^ cos Ok). 

Remark 5.3. (i) The pair (Sl,Sl) takes values in the unit circle {y^ + z"^ = 1}. 
It is elementary to see that there exists Ci > 0 such that for any k and any point 
{y, z) on the unit circle, the conditional density of (S^_,_ 2 , Sl_^_ 2 ) with respect to the 
uniform probability measure, given {(S^,S^) = {y,z)}, is bounded below by ci (note 
that the claim is about the distribution of {Sl^ 2 y ^k+ 2 )y i^k+iy ^k+i))- coupling 
argument now easily implies that the process {(S^, S^), fc > 0} is mixing and converges 
exponentially fast to a stationary distribution (which is necessarily uniform) on the unit 
circle. 

(ii) The process {{Sf), fc > 0} is a random walk. Let drC = {(x, y, 2 ;) G : y'^+z'^ < 
1, X = 0}, iVs = inf{n > 0 : > s} and let Yg G drC denote the point where the light 

ray crosses drC. It follows easily from fl5.8p that the exit time Ng goes to inhnity as 
s —?■ 00 . This and (i) easily imply that the exit distribution on drC is rotationally 
invariant. 
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Let Xk = denote a step of the random walk {S'^}. We will analyze the 

distribution of Xk. By fld.dh . 

2 cos 0fc cos sin 0fc 

(5.6) Xk — . 2 ^ — ■ 2 ^- TrX- 

sin sin 0fc + cos^ 0^ 

In order to simplify notation, we dehne 14 = sin0fc. By fl4.ip . 14 has the distribution 
1X[—1,1]. Since the distribution of 0^ is supported on [—7r/2,7r/2], cos0fc > 0 and hence 
cos0fe = — V^. For the same reason cos<I>fc > 0. This implies that 

214^1 - Vj^ cos ^k _ 214^1 - ffcos$fc 


1 — 14^ cos^ <hfc 


I 4 ^sin^$fc + 1-14 

Lemma 5.4. (a) {Xk > 0} = {I 4 > 0} and {Xk < 0} = {I 4 < 0}. 
(b) E[Xk] = 0, E[|Xfc|] = 2 - 4/77, E[Xi] = 7i/2. 

Proof, (a) This part follows from (15.Tj) . 

(b) We use (15.61) to see that, 

Ep.ii= r 


[h = sin 0] = 


J-it/2 J—k/2 sin^ 0 sin^ 9 + cos^ 9 tt 

1 8 cos 6* cos 0 sin 0 

0 Jo sin^ 0 sin^ 9 + cos^ 9 
8 cos 0 sin 0 


vr 

1 

TT 

8 

TT 

4 

TT 




Iq Jo sin^ 9 + cos^ 9 

-1 / cos6' 

/ tan ^(tan6') - d9 

Iq 2 

/•7r/2 

/ 9 COS 9d9 = 2 — 4/ TT. 


dh’^dO 


Since Xk is symmetric we have EX^ = 0. A similar calculation yields EX| = 7r/2. □ 

The following dehnitions are analogous to (13.ip . Recall that Sq = 0, S'® = 

and let Zt = S'®+ and Z7 = S'®, where T4 = infin : S'® > 0} and TP = inffn : S'® < 
-*1 -^1 

0}. It follows from Corollary 13.21 and Lemma [5.41 fbi that 

(5.8) E|Z+| = El-Zf] = 044/2 = 0F/2. 

Proposition 5.5. Fort E [0,1], 


(5.9) 


lini P ! 

-5—^00 V o “I" C/c 


< i = r(«) := 


Proof. We showed in Lemma 14.41 that f/(s, f) = P . ^ 

to the Wiener-Hopf equation. 

We have from Theorem 13.281 


E[{t{Uo + Oo) - Uo) 

E[Oo] 

< f ) is the minimal solution 


P 


4, 


Us + Os 


< f I = P I Xat. > 


s-4 


Af,-1 


t 
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where 


= / h{s-y)V+{dy), 


J —oo 
oo 


ft(s) = / P ( A'l > -—^ I 




fc=i 

‘ = PrXjVo> 




= P(t([/o + Oo)-[/o>s). 


Since s i-A P(t(t/o + Oq) — Uq > s), (13.6p implies that 

lim P (= r + Oo)-Uo>s- y)lJ^{dy) 

= j" nm + Oo) -Uo> y)dy 

_ E[{t{Uo + Oo) - f/o)+] _ E[{t{Uo + Oo) - Uo)^] 

E[^i+] E[Oo] 

□ 

Lemma 5.6. The function F : [0,1] —?■ M defined in (15.9p has the following properties. 

(a) V is a continuous, increasing and convex function. 

(b) _ 47r-3/2^^ ^ Y(^t) < t (note that 2n~^0 _ ~ o.410j. 

(c) r(0) = 0, r(l) = 1, r'(0) = 2 ti -^0 _ and r'(i) = ^E[Oo + Uo\. 

Proof, (a) Since 17o and Oo are non-negative, it is clear that F is an increasing and 
continuous function. A function of the form t i-A (at — 6)^ is a convex function for 
non-negative a and h. Since t —)■ E[(t(17o + Oo) — Uo)^] is the expected value of convex 
functions, it is convex. 


(b) By the dehnition, F(t) = ^^E[(tOo — (1 — t)Uo)^]. Since —(1 —t)Uo < 0 we have 
r(t) < j^E(tOo) = t- On the other hand, {17o = 0} = {Xi > 0}, a.s., and on that 
event we have Oq = Hence, 


r(i) - “ (1 “ ^)Oo)’^(l(!7o=0) + l(f/o>0))] 

“ “ (1 “ t)Ho)+l(f/o=o)] = ^^E[(tOo) + l(f/o=o)] 


EO, 


EOn 


(5.10) 
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It follows from the symmetry of Xi and Lemma |5.4| (b) that E[Xi i(x.>o)] = ^E|l-nii = 
1 — 2 / 71 . Corollary 13.21 and Lemma [5.41 (b) imply that E[Oo] = EZ^ = 

(5.11) E[Xil(x,>o)] ^ _ 4^-3/2_ 

^ ^ EOo 

This and (I5.10p imply part (b). 

(c) It is clear that r(0) = 0 and r(l) = 1. For the derivative at f = 0 we have: 

, , r(t)-r(o) , r(t) 

r (0) = hm -= hm - 

t^>o+ t t^o+ t 


= lim 
= lim 


t^o+ fEOo 
1 


lE[(tOo — (1 — t)f/o)’^(l(Uo=0) + l(Uo>0))] 

1 


t^o+ fEOo 
IE[Xil(Xi>o)] 
EOn 


E|(( 0 „)n, Um 


E[(tOo — (1 — t)Uo) 1([/q>o)] 


+ hm 


EOq 


-E 


Oq - - — Uq 


dUo>0) 


Note that (Oq - V^o)"^ l(Uo>o) < Oq and limi^o+ {Oq - ^Uq)^ 1(c/o>o) = 0. Hence, 
by dominated convergence and fIS.lip . 

E[Xil(Xi>o)] ^ 27r"^/^ - 47r"^/^. 


r(o) = 


EOn 


-On 


-On - Ur 


For the derivative at f = 1 we have: 

F'(l) = lim = lim ^E 

1 — f EOo 

Since t i-A y^Oq — {j^Oq — Oq)^ is an increasing function, the monotone convergence 
theorem yields 




□ 


Let Ys = (Y^, Y^, y/) G drC denote the point through which the light ray exits the 
cylinder (see Figure [9]) and let Ng = inf{n > 0 : S'* > s}. Note that 


n = 0, 


QX 

'^Ng-1 


( qy _ qy \ \ qy 


S-S- 


Ns-l 


• - ^A.-l) + ^A.-l 


It is elementary to derive the following formula from the above dehnition of Yg 

Ug 


(5.12) 


U, + 0, 


< < = {W - (1 - + (y/ - (1 - t)s%,_,f < e} 


Let D,(yQ,za) = l(y,z) : {y- {1-t)y„f + {z- {1< t‘}. Note that D,(ya,zo) 
is a disc with area The definition of Di(yo,Zo), (I5.til) and 11.5.1211 yield 

(5.13) lim P(n e D,(Sl_„ Si _,)) = r(t). 
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Figure 9. Exit point Yg is marked with (•). 


Proposition 5.7. Let Leb(y4) denote the Lebesgue measure on drC. 
(a) 


(b) 


lim 


ITTl - - -^- = C)0 

lim,^ooP(n e ] 


Leb(a.C'\A(^^,_iAk-i)) 


27rE[0( 


-E[Oo + Uo]. 


Proof, (a) We have Leb(A(*S'^ = vrt^ so part (a) follows from Lemma [5]6] 

(b) and fl5.13p . 

(b) By (EH, 

iim,^ooP(Wea,c'\A(^VD^Y.-i)) i-r(f) 

i Leb(a,C' \ Dt{Sl., S%_,)) 1 4?- 71(1 - P) 


= lim 


1 r(i)-r(«) r'(i) 


7r(l + t) 1 — t 


2 tt 


Part (b) now follows from Lemma [5.61 fcL 
Let A(0) = {{y, z) G drC : y'^ + z'^ < r^}. 

Theorem 5.8. For r G (0,1), 

(5.14) lim P(W e A(0)) < r f - t( 4^ 

s^oo \ Z / \ 2 

(5.15) lim P(n G A(0)) < a{r) := (1/2 + _ 2TT-^/^)r + 1/2 - ^ 


□ 


where V is given by (15.91) . 

Remark 5.9. The linear function a(r) takes values a(0) ~ 0.29 and a(l) = 1. 
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Proof of TheoremlEIE Since ^r-(O) C D(i+r)/ 2 iyo, zo) \ D(i-r)/ 2 {yo, z p) for any ( i/p, zp) 
on the unit circle, we obtain (15.14^ by applying (I5.13p . We derive (Ib.lbh from (I5.14p 
and Lemma 15.61 (b). □ 


5.1. Brightness singularity. We will show that the apparent brightness of the light 
arriving at the eye placed at the center of the tube opening drC goes to inhnity close 
to the center of the held of vision as the light source moves to inhnity. The precise 
statement of the result is the following. 

Let = {Sns — Sns-i)/\Sns ~ >S' 7 Vs-i| be the unit vector representing the direction of 
the light ray at the exit time. Let !B(r) = {(x, y, z) : x'^+y^+z'^ = 1, y‘^+z^ < r^,x > 0} 
denote a ball on the unit sphere and recall that Br{0) = {(y, z) G drC : y"^ + z"^ < r^}. 


Theorem 5.10. For any 0 < ri < r 2 < 1, 

hm hm (v, G ® \ ® , W G 

s^cxD 7r6^ V Vs/ Vs/ / 


r2 - ri 
2'k'^ 


The proof of the theorem will be preceded by a lemma. The lemma is an estimate 
for the Green function of the random walk Sf. The estimate is rather standard and it 
is likely to be known but we could not hnd a ready reference. Let Ms{xi,X 2 ) be the 
number oi k < Ng — 1 such that Sf — s & {xi,X 2 ). 


Lemma 5.11. For any 0 < oi < 02 < 1, 

lim —KMs{—sa 2 , —sai) = —— 

s—)-oo 7T 

Proof. Let Z1[0, 00 ) denote the space of ROLL functions equipped with the Skorokhod 
topology. Some of the functions in this space can be “killed.” We formalizing this 
idea by adding a “coffin” absorbing state to the state space and sending there all killed 
functions. We will use the convention that all functions take value 0 on the coffin state. 

Let {Wf,t > 0} be the one dimensional Brownian motion with = —s. Let 
To = inf{f > 0 : Wf = 0} and let Gs{x) denote the Green function of killed at time 
To, i.e., Gs is the function dehned by the requirement that for all —cxo < xi < a ;2 < cxo, 

PX2 pro 

/ Gs{x)dx = E 
Jxi Jo 

It is standard to show that 


—s X < —s, 

(5.16) Gs{x) = < —X —s<x<0, 

[0 x > 0. 

Recall from Lemma 15.41 (b) that EXk = 0 and EX^ = 7r/2. According to the 
Skorokhod embedding theorem (see |Obl04] ). there exist stopping times T^, k > 0, 
such that To = 0, {Tk — Tk-i,k > 1} is an i.i.d. sequence with E{Tk — Tk-i) = n/2, 
and > 1} are i.i.d. with the same distribution as {Xk,k > 1}. Let 

= inf{/c > 0 : W^{Tk) > 0} and let (xi, X 2 ) be the number of A; < — 1 such 

that W^(Tk) G (xi,X 2 ). It will suffice to prove the lemma for in place of Mg. 
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Fix an arbitrarily small £ > 0. Let Ci > 0 be so large that, for all k, a.s., 


(5.17) 


E 


'n- 




Suppose that 0 < oi < 02 < 1 and s is so large that {—sa 2 — ci, —sai + ci) C (—s, 0). 
Then fl5.16p and (15.1711 show that 




(7r/2 - e)'EMY (-soa, -soi) = E ^ {lw’>{n_^)&{-sa2-sa^){'^/'^ - £)) 


k=l 


Ee(i< 


^k<N'^ -*-IT'’(Tfc_i)e(—sa2,—sai) 

k=l 

00 

< ^ ^ E I \h<N'w'\. 


(7r/2 - £)) 


k=l 


k<N^ -*-IF'’(Tfc_i)e(—sa2,—sai) 
Tk 

Tk-i 


X E / l(^w‘‘{Tk-i)-ci,W‘>{Tk-i)+ci){W^)dt I 


rn 


j lk<N^lw’>{Tk-i)e{-sa2-sai) l(lT'>(rfc_i)-ci,IT'>(Ti,_i)+ci) (hF/)dt 


k=l 




— ^ I 1(—sa 2 —Cl,—sai+ci) (hF,f )dt 


r'^o 




<E / 1(_SQ,2-Ci,-Sai+Ci) (hFt )(it + E / 1(_SQ,2-Ci,-Sai+Ci) (hFt )(it 


o—sai+ci 


< 


Gs{x)dx + E {Tjsfw — Tjsfw_i) 


' —sa2—ci 


= ^{{-sa 2 + cif - {-sai - cif) + 7r/2. 


It follows that 


1 1 ~ ®l) 

limsup —EMs(-sa 2 , -sai) < - —-■ 

s—>00 S j Z E 


Since £ > 0 is arbitrarily small, 


(5.18) 


limsup —EMs(—soa, —sai) < 


(4 - aj 

TT 


Let 

8^7172 

The processes {Vf,t > 0} converge to {Wl,t > 0} in distribution in the Skorokhod 
topology when s —)■ cx). 
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For a, b and 5 > 0 such that a+S < b—6, choose some continuous function A(a, b, S, t) : 
R —)■ [0,1] such that 


A(a, b, 6, t) 


{ 0, if t < a, 

1, iia + 6<t<b — 6, 

0, if t > b. 


Fix any u < oo. The functional / /“ A(a, b, S, f{t))dt is bounded and continuous on 

iA[0, oo) in the Skorokhod topology. It follows that 


lim inf E 

S^OO 



Ms{—sa2, -soi) 


> lim inf E [ A f— 02 /\/ 6, dt 

s^oo Jo ^ ' 

= E f X dt 


> 


™—ai ^ 'K l‘l—6 


Gs{x)dx 


— “ 2 /^/ 7r/2+<5 

(^-a2/V^/2 + 5^ - {-ai\/n/2 - <5^ ^ . 


Since 5 > 0 can be arbitrarily small, 

lim inf E Ms(— 502 , —soi)^ > — - 

s—>-oo y TT 

This and flb.lSh prove the lemma. □ 

Proof of Theorem \5.1(A Recall that 5^(0) = {{y,z) G drC ■. < r^}. If s is 

large, 5 > 0 is small and the light ray leaves a point {—sa, 0,1) in the random direction 
determined by fl5.3l) - fl5.5l) then the probability that this ray will exit the tube through 
^ 5 ( 0 ) is 


(5.19) 


7r(5^ 


1 1 

27rsa 2(sa)^ 


+ o{6‘^/s^) 




+ o((5^/s^). 


The factors on the left hand side represent the following quantities. The area of Bs{0) 
is equal to n6‘^. The hitting density is the product of two factors, corresponding to $ 
and 0 (see the beginning of Section [5] for the dehnitions). The factor representing the 
density of $ is (the reciprocal of the radius of the circle centered at the starting 
point and passing through the center of drC, up to a term of lower order). The factor 
representing the density of 0 is 2 {L)^ because of fl2.ip and scaling by the radius sa, just 
like in the case of the density of <F; once again, the terms of lower order are ignored. 

For hxed ri and r 2 , large s, and small 5 > 0, a light ray arriving at Bs{0) in the 
direction G R (y) \ ® (^) niust have have left the surface of the tube at a point 
with the x-coordinate in the range (—s/ri + o(s) + 0{6s), —s/r 2 + o(s) + 0((5s)). 
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We define a measure by Ms(y4) = E ( J2k=o ^ ) for every Borel subset A of 

Lemma 15.111 can be rephrased as 


lim — 

s—)-oo 


,{{—sa2, —soi)) — 


gj - 

71 


A formal calculation based on this formula yields for small e > 0, 

2aes^ 


(5.20) 


,(-s(a + e), -sa) 


vr 


It is routine, using techniques from the proofs of Lemmas 13.31 and 13.211 to provide a 
rigorous argument based on fl5.19l) and fl5.20p . showing that for any 0 < ri < r 2 < 1, 


lim lim G ® ) \ ® f—) , W G ^^(O) 

s^-oo 5^0 V Vs/'-'' 


= lim lim 


c--l/r 2 +o{s)/s+ 0 {ds)/s 

r (52 , , 1 

1—1 /ri +o(s) /s+0{Ss) / s 

3 + o{5^/s^) 
4s3a3 



da = 


r2 - Ti 
271^ 


M(s da) 


□ 


5.2. Discussion of the results. The behavior of the light reflection process in the 
three dimensional tube is much different from that in the two dimensional case. The 
most notable difference is that the overshoot Os and undershoot Ug (in the x-direction) 
converge to a non-trivial distribution (instead of going to infinity as in fl4.5p ). The 
reason is that the ladder variable has finite expectation (see (I5.8h ). unlike in the 
two dimensional case. This fact and the Wiener-Hopf equation can be used to show 
existence of the limiting distributions for many quantities of interest. Unfortunately 
most of the formulas that can be obtained in this way are abstract integrals that cannot 
be easily interpreted. 

Theorem 15.101 savs that small annuli at the center of the held of vision, with the area 
of magnitude 1/s^, receive about 1/s units of light. Hence, the apparent brightness is 
about s at the distance 1/s from the center, if the light source is s units away. This 
means that the surface of the tube does not appear to the eye to be Lambertian, i.e., 
the surface does not have uniform apparent brightness. This can be explained by the 
fact that not all parts of the surface of the tube receive the same amount of light. 
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